Soliton solutions of driven non-linear and higher 
order non-linear Schrodinger equations 



Vivek M Vyas^, T Soloman Raju ^, C Nagaraja Kumar'^ and 
Prasanta K Panigrahi^ f 

^ Department of Physics, M S University of Baroda, Vadodara, 390 002, India 
^ Physics Group, BITS-Pilani, Goa Gampus, Zuari Nagar, Goa, 403 726, India 

Department of Physics, Panjab University, Chandigarh, 160 014, India 
''Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India 

Abstract. We analyse the structure of the exact, dark and bright soliton solutions of 
the driven non-linear Schrodinger equation. It is found that, a wide class of solutions 
of the higher order non-linear Schrodinger equation with a source can also be obtained 
through the above procedure. Distinct parameter ranges, allowing the existence of 
these solutions, phase locked with their respective sources, are delineated. Conditions 
for obtaining non-propagating solutions are found to be quite different for both the 
equations. A special case, where the scale of the soliton emerges as a free parameter, 
is obtained and the condition under which solitons can develop singularity is pointed 
out. We also study the highly restrictive structure of the localised solutions, when the 
phase and amplitude get coupled. 



PACS numbers: 42.81.Dp, 47.20.Ky, 42.65.Tg, 05.45.Yv 
1. Introduction 

The externally driven, non-linear Schrodinger equation (NLSE) with a source has been 
investigated in the context of a variety of physical processes. It arises in the problem 
of Josephson junction [1], charge density waves [2], twin-core optical fibres [3, 4, 5, 6], 
plasma driven by rf fields [8] and a number of other problems [9]. As compared to NLSE, 
which is an integrable system [10], not much is known about the exact solutions of this 
equation. Perturbative solutions around the stable soliton solutions of NLSE with a 
source have been studied earlier. Analysis around constant background and numerical 
investigations [11, 12, 13] have revealed the phenomenon of auto-resonance [14, 15] as a 
key characteristic of this system, where a continuous phase locking between the solutions 
of NLSE and the driven field is observed. 

In a recent work [16] , some of the present authors have devised a procedure based on 
fractional linear transformation for obtaining exact solutions of this dynamical system. 
One obtains both localised and oscillatory solutions. Apart from these regular solutions, 
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under certain constraints singular solutions have also been found, implying extreme 
increase in the field intensity due to self-focussing. This approach is non-perturbative in 
the sense that, the obtained exact solutions are necessarily of rational type, with both 
numerator and denominator containing terms quadratic in elliptic functions. However, 
the exact parameter ranges in which the general solution exhibits bright and dark nature 
has not been investigated. Considering, the importance of the localised solutions of this 
physically important dynamical system, the above point needs a systematic study. 

The goal of the present paper is to analyse in detail the structure of the most 
general, locahsed bright and dark sohtons of the NLSE with a source. The parametric 
restrictions, under which singular structures can form, are obtained. Unlike NLSE, 
it is observed that dark and bright solitons, depend both on coupling and source 
strengths. For example, bright solitons can also form in the repulsive regime, if the 
source strength has positive value. Similarly dark solitons can form in the attractive 
regime. Conditions which give non-propagating solutions, for driven NLSE, are studied. 
The possibility, where the phase and amphtude can get related is also investigated. The 
highly restrictive nature of the resulting dynamics is pointed out. In the process, we also 
find a large class of solutions to higher order nonlinear Schrodinger equation (HNLSE) 
with a source, by connecting the solution space of these two equations. HNLSE, which 
is a generalisation of NLSE, has been proposed by Kodama [17] and Kodama and 
Hasegawa [18] to describe the propagation of short duration (femtosecond) pulses in 
optical fibres. The fact that NLSE with a source describes twin-core fibres, under 
certain conditions [5, 6, 7], motivates us to find exact solution of HNLSE with a source. 
A possible scenario leading to driven HNLSE is pulse propagation through asymmetric 
twin-core fibre involving a femtosecond pulse in first core, having anomalous dispersion, 
and a nanosecond pulse propagating in second core in the normal dispersion regime. 
Analysis of the solution space of the driven HNLSE shows many interesting features. 
Non-propagating solutions of driven HNLSE show a much richer structure as compared 
to NLSE, because of the presence of the higher order terms. 

The paper is organised as follows. In the following section, we elaborate briefiy on 
the method of the fractional linear transform for obtaining the solutions of the NLSE 
with a source. In this section, it is shown that for a large class of solutions, driven HNLSE 
can also be connected with driven NLSE. Hence, the obtained solutions are also solutions 
of HNLSE, albeit with different relations among the parameters. The third section is 
devoted to the investigation of appropriate parameter regimes, characterising bright, 
dark and singular solutions, and their analysis. In the fourth section, the structure of 
the solutions space is studied, when the phase and amplitudes are coupled. We then 
conclude after pointing out the open problems and future directions of work. 

2. Analysis of NLSE phase locked with source 
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and its relation to HNLSE 

The equation which we intent to solve is driven NLSE, which is space time dependent, 
and phase locked with the source: 



86 d"^ 

H 

dt dx^ 



+ ^ + 9 \ cf> cf> + f^<P ^ '^e'^"''-^'^ , (1) 



where g, /i, k, u> and k are real constants. As we will see later, it can be connected with 
HNLSE, for a wide class of solutions. We consider the ansatz traveUing wave solution 
in the form, 

where ^ = a{x — vt). Separating the real and imaginary parts of equation (1), one 
obtains, v — 2k, from the imaginary part, indicating that in the present case, wave 
velocity v is controlled by k. The real part yields, 

a'^a" + ga"^ + ea = k, (2) 

where e = u — k^ + fi, and prime indicates differentiation with respect to ^. As has been 
observed earlier [16], this equation can be connected to the equation /" + a/ + bf^ — 
through the following fractional linear transformation (FT): 

"«> = 1 + D/K|m)= 

where A,B and D are real constants, and is a Jacobi elliptic function, with the 

modulus parameter m. We consider the case where, f{^\m) = cn{^\m), other cases can 
be similarly studied. Since the goal is to study the localised solutions, we consider the 
case with modulus parameter m = 1, which reduces cn{^) to sec hyperbolic of C,- It 
is worth pointing out that, other solutions involving sn(^) and dn{^) naturally emerge 
from the above solution, since the transform involves square of the cn{C,) function. It 
should also be mentioned that, in the above fractional transform the second power of 
the cnoidal functions emerge, as the highest power. 

We can see that equation (3) connects a{^) to Jacobi elliptic equation, provided 
AD^B, and the following conditions are satisfied for the localised solutions: 

Ae + gA^ - 0, (4) 

2eAD + eB + W{B - AD) + 3gA^B - 3kD = 0, (5) 
AeD^ + 2eBD + W{AD - B)D 

+6a^(AD - B) + ^gAB"^ - SkD'^ = 0, (6) 

eBD^ + 2a^(B-AD)D + gB^-KD^^0. (7) 

Equation (4) in A does not involve B and D, which is first solved to get the real A. 
Thus A is determined in terms of e, k, and g. From equation (5), we determine D in 
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terms of S as, D = FB, where, T = 4a2^A+3K-2€A • substituting this in equation (6), 

B is found as -B = 3g^+^gr2+2er+4a2r(!4r-i)-3K;r2 • equation (7), we obtain a cubic 

equation in (3 = a^: 

Pi/?^ + gi/3' + ri/3 + ci = 0, (8) 

where pi = 64:{A^g - Ae - k), qi = {ASA^g'^ + BAA^ge + 16 Ae'^ - ASA^gK - IGe/i), 
n = (UA'^g^ + 36A^g^e + 20A^ge'^ - AAe^ - GOA^g'^n - nA^geK + Ae^K + Ad>AgK'^) and 
ci = (3AVe - "iA^g^e^ - 7A^ge^ - Ae^ - 18AV« - IS^Ve^ + I'^A^ge'^i^i + + 
^^g^K^ — IbAgen^ + 9gi<^). It can be straightforwardly seen that pi in equation (8) 
is the consistency condition (4) and hence is identically zero. Therefore, the width 
parameter [3 is the solution of a quadratic equation. Thus for any given values of g, e 
and K,, wc can find the values of A, B, D and a. 

Before continuing the analysis of the localised solutions, let us turn to driven 
HNLSE: 



dx^ dx dx 



where g, n, e. A, B, C, lu and k are real constants, and k is the source strength. 

The different nonlinear terms describe various interactions which affect the 
propagation of femtosecond pulses in optical fibres [19]. The term proportional to A 
results from the inclusion of the effect of third order dispersion and the term proportional 
to B comes from the first derivative of the slowly varying part of nonlinear polarization, 
which is responsible for self-steepening and shock formation. The delayed Raman 
response for self-frequency shift accounts for the term proportional to C. 

In last many years lot of activity has gone into the study of the solvability of 
the above equation without a source. Integrability studies were performed using inverse 
scattering transform method and Hirota method and a few cases of integrability has been 
identified:(i) Sasa-Satsuma equation {A:B:C — 1:6:3) [20], (n) Hirota equation {A:B:C 
— 1:6:0) [21], (in) derivative NLSE, of type I and type II [22]. Several exact solutions of 
solitary wave (bright sohton) and kink type (dark soliton) have been obtained [23, 24]. 
Phase modulated solutions to HNLSE have also been obtained [25], where the generic 
form of the solution reads 'i/j{x,t) = /(a;, t)e**--^'-^'*'*~^*^; here x is some function of x and 
t. Substituting the ansatz solution of the form, 

^{x,t)^p{Oe'^''^-^'l (10) 

in equation (9) and then separating out the imaginary and real parts, we get the following 
two equations, 

a^eAp" + {2ak - 3aek'^A - av)p' + {aeB + 2aeC)p^p' = 0, (11) 
a^p" + gp^ + ep^ K. (12) 



Here the parameters are defined as, 

^2 = - 3eAk), 



(13) 
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g^{g-eBk), (14) 

and 

e = {u - P + fi + eAk^). (15) 

It is interesting to note that, the above real part of HNLSE is similar to the real part 
of driven NLSE. 

Differentiating equation (12) once and comparing it with equation (11), we see that they 
are consistent only if the parameters satisfy the relations: 

5 + 2(7=^ (16) 

and 

v^2k-3eAe-^S^. (17) 

Hence, for this class of solutions the problem of solving equation (9) is mapped to 
the problem of solving equation (12), which is nothing but the real part of the driven 
NLSE. This proves our assertion that phase locked solutions of driven NLSE are mapped 
to the phase locked solutions of driven HNLSE, in certain parameter ranges. One can 
see from the above relation that the velocity in this case has a quadratic relation with 
k, a feature very different from the driven NLSE. 

3. Analysis of the localised solutions 

We now proceed to analyse carefully the localised solutions of equation (2). It is clear 
that a localised solution is bound to have at least one extremum in its profile. This 
implies that the first derivative of a{^) must vanish at the extremum: 

2(B-AD)m 

(1 + omy 

Since AD^B, either / or /' or both must be zero. In our case /(^) = sech{^), whose 
first derivative vanishes only at origin. This means we have an extremum at origin. The 
second derivative at origin is, 

2(AD-B) 

which resolves the maximum and minimum. It should be noted that a" is singular for 
D — —1. For the non-singular case, we see that there is a clear distinction of two 
regimes of solutions: One for which AD > B, where a" is positive; this corresponds to 
minimum. In the second case, where AD < B is negative, we have a maximum. The 
latter, corresponds to a bright soliton, whereas the former corresponds to a dark soliton 
in the propagating media. This clearly suggests, that both types of solitons exist in this 
dynamical system. 
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Figure 1. Contour plot showing the allowed background field amplitude value A, for 
which the localised solitary wave solution exist, here e = 0.1. 



In these rational solutions, parameter A decides the strength of the background, in 
which these solutions propagate. It is interesting to note that for the locahsed solution, 
A — Q is not permitted, since this leads to the absence of the source. Considering a 
to be positive scmi-dcfinite one finds further constraints: (1)^4 should be positive since 
negative background is not meaningful; {}i)A > — (iii)-D should be greater than — 1 
for non-singularity. The parameters satisfying these conditions are taken to be physically 
meaningful. The parameter conditions leading to singular solutions are dealt separately. 
Figure (1) illustrates various ranges of the values of the strength of the nonhnearity g, 
and strength of the source re, for solitary wave solutions. Figure (2) shows the solutions 
for some mentioned values of the parameters. 

In case of driven HNLSE, we see that v — 0, when the equation, 

Se'^A^k^ - SeAk'^ + 2k- €A{w + n) ^ 0, (20) 

is satisfied. In general, the above cubic equation can have one real, non-zero value of 
k, for appropriate equation parameters. Hence, we see that driven HNLSE has non- 
propagating solutions even when k^^O, which is different from driven NLSE, since in the 
latter case v — 2k. 

Considering the case ^ = 0, we see that, equation (12) and equation (2) have same 
coefficients and hence the same solutions. However, unhke the driven NLSE for this 
case one can have a propagating envelope with velocity 

V = -eeA. (21) 

Static solutions are obtained, when any one of A, e and e vanish. It is interesting to 
note that when e vanishes, as we shall show later, a becomes a free parameter, hence we 
get a static sohton of arbitrary size. The case when e vanishes, driven HNLSE simply 
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Figure 2. Density profiles of some solitons of driven NLSE. 

i) Dark soliton (dashed), with g = 0.7, k = —0.7 and e = —1.4; 

ii) bright soliton (solid), with g = —0.5, k = 0.1 and e = —0.4; 

iii) bright soliton (dash-dotted), with g = —0.4, k = 0.1 and e = —0.4. 



reduces to driven NLSE, for which k = imphes v = 0. Considering the case when 
the third order dispersion parameter A vanishes, one finds that for the non trivial case 
B -.C = -2:1. 

We now consider the case when e = for which uj — k'^ — fj,; in this case equations 
(4) to (7) yield, 

qi{a^)+ria'^ + c = 0. (22) 

Here, gi = ASA^g'^ - ASA^gK, ri = 12AV - G0A'^g2K + ASAgK^ and c = -ISA^g^n + 
dA^g'^K^ + 9gK^. Very interestingly, all of these coefficients vanish in the view of equation 
(4), leaving a as a free parameter. So, the width of the soliton, in this case, is 
independent of the parameter values, which means that the solitons can have arbitrary 
size for the given values of g and k. Similar is the case for driven HNLSE, when e 
vanishes from equation (12). 

As noted earlier, equation (2) has singular solutions, for D — —1. One can show 
that they satisfy the relation 

g{6a'^ - 3e)A^ + e{2a'^ + e) + 9gKA - = 0. (23) 

These solutions possess one singularity in their profile, which physically corresponds 
to very large field intensity due to self focussing [26, 27]. When D < —1, singular solitons 
exists, with singularites existing at two different locations, between which a becomes 
negative. 
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4. NLSE with a generalized source 

Instead of the source term of the type we have considered previously, one can have a 
more general one where, 



^^ + ?T + ^lV'rV' + /^V^ = Ke^(>^«)-'^*) . (24) 
ot ox^ 

Here xH) is some function. We consider the ansatz, 

i;{x,t)^p{i,t)e'^^^^--'\ (25) 

where ^ = a{x — vt). Equation (25) when substituted in equation (24) gives a complex 
equation in p and x- Equating the imaginary part to zero, one gets, 

^' = ^ + ^- ^26) 
la p'^ 

One can clearly see that the above equation suggests phase- amplitude coupling for c^Q. 
We see from relation (26) an interesting case of the phase singularity of the soliton 
arising through phase-amphtTidc couphng, when p vanishes. 

Substituting this relation in the real part of equation (24), we get a non linear differential 
equation of the form, 

o?p" + ^ ep = K + — - , (27) 

where e = {u -\- p -\- '^). The above equation can be connected with the equation 
governing Jacobi elliptic functions via a fractional transformation as, 

_ A + Bf{i\mf 

- 1 + D/KW" '''' 
where 5 — 2\s again the maximum allowed value. This gives the consistency conditions: 

gA"^ + eA" - {2a\AD - B){1 - m) + k)A^ = a^c" (29) 

6A^Bg + 4eA^B + 2eA^D - ZkA^B - ZkA^D + Qa'^iB - AD){1 - m)A^B + 
W{AD -B){1- m)A^D + 4a^(2m - 1)A%B - AD) = Qa^c^D (30) 

6a^(S - AD){1 - m)AB^ + 18a^{AD -B)(l- m)A^BD + 
l2a^{B - AD){2m - l)A^B + W{AD - B){2m - 1)A^D + 
da^iAD - B)mA^ + IbgA^B'^ + GeA^B^ + %eA^BD + eA^D'^ 
-3kAB2 - QkA^BD - SkA^D^ = IBa^c^D^ (31) 
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-kB^ - 9kAB^D - QkA^BD^ - kA^D^ + 

18a'^AB^D{AD - B){1 - m) - 2a^B\AD - B){1 - m) - 

12a^AB^{AD - B){2m - 1) - 2a^mA^D{AD - B) 

+l2a^A^BD{AD - B){2m - 1) = 2Qa^c^D^ (32) 



IbgA^B"^ + eB^ + %eAB^D + 6eA^5^D^ 

-SkB^D - 9kAB^D^ - SkA^BD^ + 6a'^B^D{AD - B){1 - m) 

-Ao^-^B^AD - B){2rn - 1) + 12a^AB'^D(AD - B)(2rn - 1) + 
18a^mAB'^{AD - B) - Qa^mA^BD{AD - B) = iWc^D^ (33) 



QgAB^ + 2eS^D + AeAB^D'^ - SkAB^D^ 
-SkB^D^ + Aa^iB - AD){1 - 2m)B^D 

+6a'^m(B - AD)AB'^D + Q>a^m{AD - B)B^ ^ da^c^D^ (34) 

2a\B - AD)mB^D + gB^" + eS^D^ - kB^D^ = a'^c^D^ (35) 

Since, there in total are seven simultaneous relations to fix three independent 
parameters, we can see that these are constrained solutions. One can also see that the 
phase-amplitude coupling imposes constraints on the solution space, without changing 
the profile of the solutions. Prom the above relations, one can see that 5 = needs 
c = or D = 0. The former case has already been studied here, whereas the latter case 
results in a constant background solution. Here, A^Q requires ct^O, hence the solutions 
always exist in a constant background. This should be contrasted with the case when 
c = 0; for which cnoidal wave solutions are possible with A = 0. Similarly, in m = 1 
case localised solutions with B — {) are allowed in the previous case, which as seen above 
are not found in present case. Hence, only rational solutions are possible. Moreover, 
we see that (29) is a sixth order polynomial, in contrast to the previous relations, this 
does not have analytically tractable roots leaving numerical analysis as the only tool to 
analyse the structure of solution space. 

In conclusion, a number of interesting features have emerged from analysis of the 
exact solutions of the driven NLSE and HNLSE. Dark and bright solitons can exist 
in attractive and repulsive non-linear regimes, a feature very different for NLSE; the 
presence of the external source makes this possible. In a wide range of parameters, the 
driven HNLSE and driven NLSE have similar solutions, albeit with different sizes and 
velocities. Static solitons are found for both the equations. In particular, static solitons 
are found to have interesting properties like arbitrary scale, under certain parametric 
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restrictions. They differ significantly for driven NLSE and HNLSE cases. In certain 
specific parameter regimes, solitons of arbitrary size are found. Singular solutions are 
also found for both the equations. The investigation of the situation, where the phase 
is allowed to depend upon the intensity, revealed that the corresponding solutions are 
highly constrained. The study of solitary waves having complex envelope, analogous 
to Bloch solitons in condensed matter physics, is worth investigating in the present 
scenario. Application of the fractional linear transformation technique employed here 
to other nonlinear equations is also of deep interest. Investigations along these lines are 
currently in progress. 
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